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Abstract 



We interpret the phase structure of hadronic matter in terms of the basic dynamical and 
geometrical features of hadrons. Increasing the density of constituents of finite spatial 
p , extension, by increasing the temperature T or the baryochemical potential fx, eventually 

X/, "fills the box" and eliminates the physical vacuum. We determine the corresponding 

D I transition as function of T and /i through percolation theory. At low baryon density, this 

'^ ■ means a fusion of overlapping mesonic bags to one large bag, while at high baryon density, 

hard core repulsion restricts the spatial mobility of baryons. As a consequence, there are 
^ . two distinct limiting regimes for hadronic matter. We compare our results to those from 

0\ ! effective chiral model studies. 

1^ : 1 Introduction 

o. 

OO , In recent years, the general features of the phase diagram of strongly interacting matter 

have become increasingly well established [1]. Lattice QCD studies at finite temperature 
and now also for some range of finite baryon density [2-5], combined with chiral synime- 

/\ • try restoration arguments [6-12], lead for physical quark masses to the phase structure 

j^ ■ illustrated in Fig. [1] as function of temperature T and baryochemical potential ft. Follow- 

ing a region of non-singular but rapid cross-over of thermodynamic observables around a 
quasi-critical temperature of 170 - 190 MeV, increasing fi leads to a critical point, beyond 
which the system shows a first order transition from confined to deconfined matter. 

The aim of the present work is to deduce such behaviour from the basic dynamical and ge- 
ometrical features of hadronic matter. Our question thus is the following: what conceptual 
aspects of hadronic interactions lead to the observed behaviour, and in particular, what 
features in hadronic dynamics result in the observed changes of the transition structure 
as function of T and /i? 

At low baryon density, the constituents of hadronic matter are mostly mesons, and the 
dominant interaction is resonance formation; with increasing temperature, different reso- 
nance species of increasing mass are formed, leading to a gas of ever increasing degrees of 
freedom. They are all of a typical hadronic size (with a radius Rh — 1 fm) and can over- 
lap or interpenetrate each other. For /x — 0, the contribution of baryons/antibaryons and 
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Figure 1: Phase structure of QCD matter in the T — fi plane 

baryonic resonances is relatively small, but with increasing baryon density, they form an 
ever larger section of the species present in the matter, and beyond some baryon density, 
they become the dominant constituents. Finally, at vanishing temperature, the medium 
consists essentially of nucleons. 

At high baryon density, the dominant interaction is non-resonant. Nuclear forces are short- 
range and strongly attractive at distances of about 1 fm; but for distances around 0.5 fm, 
they become strongly repulsive. The former is what makes nuclei, the latter (together with 
Coulomb and Fermi repulsion) prevents them from collapsing. The repulsion between a 
proton and a neutron shows the purely baryonic "hard-core" effect and is connected 
neither to Coulomb repulsion nor to Pauli blocking of nucleons. As a consequence, the 
volumes of nuclei grow linearly with the sum of its protons and neutrons. With increasing 
baryon density, the mobility of baryons in the medium becomes strongly restricted by the 
presence of other baryons, leading to a "jammed" state, as shown in Fig. [2] [13]. 
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Figure 2: Hard sphere states: full mobility (a), "jammed" (b) 



Increasing the density of constituents by increasing temperature, baryon density or both, 
leads to clustering of hadrons of spatial extension, and these clusters will eventually span 
the entire available volume. This onset of connectivity and the associated geometric 
critical behaviour is the central topic of percolation theory, where the relevant thresholds 
have been determined for permeable spheres which can overlap (mesons) as well as for 
those which have a impenetrable hard core (baryons). In particular, one finds in both 
cases a geometric transition from a state in which the vacuum forms a finite part of the 
system to one in which only isolated bubbles of vacuum remain in a world of fully or 



partially overlapping hadrons. We shall consider this point of "disappearance" of a large- 
scale vacuum as the end of hadronic matter and calculate the corresponding limiting curve 
in the T — fi plane. 

Our question therefore addresses two distinct situations. We have to consider the per- 
colation of mesons, allowing full overlap, and then that of baryons with a hard core. In 
the next section, we shall first recall the salient features of percolation theory for the two 
cases. Following this, we will use a very simplistic toy model to illustrate the underlying 
concepts. In section 4, we then turn to the realistic case of a hadronic resonance gas with 
baryon number and strangeness conservation. In the last section, we consider the nature 
of the transition and compare our conclusions to the results obtained in chiral model 
studies. 



2 Percolation 

Here we briefly recall the essential results for the percolation of spheres, in three space 
dimensions, for the case of arbitrary overlap [14] and for spheres with an impenetrable 
hard core, allowing only partial overlap [15]. Consider A^ spheres of radius Rq and hence 
volume Vo = (47r/3)i?Q in a "box" of size V, with V ^ Vq. Percolation is said to occur 
when a connected and hence at least partially overlapping set of spheres spans the volume, 
or in other words, when the volume occupied by the largest cluster of connected spheres 
reaches a finite fraction of V. 

We address first the case of permeable spheres, i.e., with arbitrary overlap. In this case, 
when the density n = N/V increases, clusters of overlapping spheres form, and for 

^m = -rT- (1) 

with rjm — 0.35, the largest cluster first spans the system, i.e., percolation occurs [14]. At 
this point, the fraction 

(j)^ = l- e"''™ ~ 0.30 (2) 

of space is occupied by spheres; the complementary 70% remain empty space ("vacuum"), 
which also still spans V. The largest cluster therefore has a density of about 1.2/Vo. In 
other words, at the percolation point, a randomly created state of a system of N spheres 
is very inhomogeneous, consisting mainly of one dense cluster and much empty space. A 
random distribution of extended constituents thus leads to something like a "geometric 
attraction" between pointlike constituents, with clustering as the result. We shall see 
later that in the percolation of hard-core spheres, this attraction is competing with an 
intrinsic repulsion. 

For n < rim, the clusters of overlapping spheres form only isolated bubbles in V. From eq. 
([1]) we see that at the percolation point, the total volume of all spheres adds up to 35% 
of the total volume; this is covered by the spheres to only 30%, indicating 5% "overlap". 
After a further increase of density, a second percolation point is reached at 

r^v = TT, (3) 

^0 



where r]^ ~ 1.22. Now 

(j)^ = l- e-"^ ^ 0.70 (4) 

is the fraction of space covered by spheres, with only 0.30 remaining empty. Here the sum 
of the sphere volumes is with 1.22 V considerably larger than the total volume, due to 
increased overlap. For n < n„, the vacuum percolates, above n„, only isolated bubbles of 
vacuum remain. In other words, for n > n„, the vacuum has "disappeared" as a large-scale 
entity. 

The existence of two percolation thresholds, one for the formation of the first spanning 
cluster of spheres and another for the disappearence of a spanning vacuum "cluster" , is a 
general feature of percolation theory in three or more dimensions. 

We now turn to the percolation of spheres (again of radius Rq) having an impenetrable 
spherical hard core [15], which we assume to have the radius Re = Ro/2. Each sphere now 
defines a volume Vq = (47r/3)i?o which is not accessible to the center of any other sphere. 
The spheres can partially overlap: the distance between their centers has to remain greater 
than Rq = 2 Re. With increasing density, we now have again two percolation thresholds. 
At 

Vm /p-N 

IT-m — 'TT' V / 

the spheres form a spanning cluster, and at 

nv = Tr (6) 

the vacuum last spans the volume. Numerical studies [15] show that f]m — 0.34, practically 
the same as found above for permeable spheres of the same size. The percolation threshold 
for spheres with a hard core is on the dilute side thus not much affected by the presence 
of the hard core, and we have a similar geometric attraction. The vacuum percolation 
threshold, however, now is given by fj^ ~ 2.0, in contrast to rj^ ~ 1.24 for permeable 
spheres. In other words, the disapperance of the vacuum requires for hard core spheres a 
higher density than needed for permeable spheres. This is directly related to the hard core 
repulsion, which tends to move the spheres apart and thus at high density counteracts 
tight clustering. 

3 A Toy Model 

In our simplified toy model, we want to compare the percolation behaviour of an ideal 
gas of massless pions of radius Rh to that of an ideal gas of massive nucleons of the same 
size, but having a hard core of a (smaller) radius Re = Rhf^- The density of pions is 
specified by the temperature T of the medium, that of the nucleons by the temperature 
T and the baryochemical potential n, fixing the overall baryon number (nucleons minus 
anti- nucleons). 

Consider first the system of pions. With three charge states, the density is given by 

nAT) = 3 ^ Tl (7) 



With increasing temperature, they will overlap and eventually fill the given volume with 

one big connected bag. When only isolated vacuum bubbles remain, we assume the system 

to have reached the limit of pionic matter. In the previous section it was shown that this 

occurs for 

1 22 
Uf = —— ~ 0.57 im-\ (8) 

where Vh = (47r/3)i?| and we have used Rh = 0.8 fm. Solving rihiT) = Uf yields 

96 

T^ ~ -— ~ 240 MeV (9) 

Rh 

as the limiting temperature obtained through pion fusion. This value will drop slightly 
when we include nucleons and antinucleons, and it will decrease considerably when reso- 
nances are brought in, as we shall see in the next section. 

The density of pointlike nucleons of mass M is at T = (when there are no anti-nucleons) 
determined in terms of the baryochemical potential /x, 

n,(/x,T = 0) = -^(/.2-M2p. (10) 

In the presence of a hard core of volume Vc = Vh/8, the density has to be calculated 
taking into account the reduction of the available volume. We approximate this by a van 
der Waals approach [16], taking^ 

neifi, T = = -— — 11 

l + nb{T,fi)Ve 

for the density of extended nucleons, where K — 2Vc denotes the excluded volume at 
random dense packing of hard spheres [18]. With increasing nucleon density, the box 
becomes more and more filled, and we saw above that the empty vacuum disappears for 

2 

n^ ~ — - 0.93 fm-^ (12) 

Vh 

taking the same radius for nucleons as for pions; this corresponds to about 5.5 times 
standard nuclear density po — 0.17 fm~'^. Solving ^^(T = 0,/i) = n^ yields 

/i^ ^ 1.12 GeV (13) 

for the limiting baryochemical potential at T = 0. 

In the region of low to intermediate /i in the T — ^ diagram, we approximate the density 
of pointlike nucleons by the Boltzmann limit 



.,.).f?(f)%«M/r)e-.?(^)" 



n,ifi,T) ^ ^ — K,{M/T)e^l' ^ — ^- e^^-^'^" , (14) 



^A thcrmodynamically consistent implementation of hard core repulsion requires in addition a shift 
of /i [17], which for simplicity is neglected here. 



where we have also used the large M/T form of the Bessel function K2{M/T). At /i = 0, 

we can use this to add nucleons and antinucleons to the pions considered so far for filling 

the box. From 

1.22 
n^ + 2nB = ^— (15) 

we then find a limiting temperature of hadronic matter 

T,,~230MeV (16) 



slightly lower than for pions alone. The small change is due to the fact that nucleons 
amount to only about 6% of all hadrons at fi = 0. If we would there want to fill the box 
with only nucleons and antinucleons, we would have to go to the much higher temperature 
of about 430 MeV. 



Using the approximation ([HD, we can ask for the value of fx at which the nucleon per- 
colation curve crosses the hadronic limit T^. Solving ub^T = T^, fi) = 2/Vh for /x, we 
find /i ^ 0.8 GeV for the crossing point. The comparitive behaviour of the two curves is 
illustrated in Fig. [31 It is evident that at low nucleon density, pion percolation limits the 
hadronic matter density, while at high nucleon density, the percolation of hard spheres 
and the resulting jamming provide the limit. It should be emphasized, however, that the 
separate calculation of a (permeable) pion and a (hard-core) nucleon curve is obviously 
not the final solution, since a dilute admixture of nucleons and antinucleons contributes in 
the low density regime, just a pions will help to reach percolation in the nucleonic region. 
We have followed an additive picture here, since so far (to our knowledge), continuum 
percolation studies have been performed either for permeable spheres or for spheres with a 
hard core. Evidently what is needed here is a study allowing both types of constituents in 
a degree of mixture specified by /i, and this will lead to modifications in the intermediate 
region, as schematically included in Fig. [31 
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Figure 3: Limit from massless pion percolation and from that of massive hard core nucle- 
ons; the curve labelled "mixture" indicates the result expected of a combined percolation 
study. 

We note here briefly that we can determine a confinement-deconfinement limit also by 
comparing the hadronic pressure to that of a quark-gluon plasma [16,19]. Such a compar- 
ison leads by construction to a first-order transition in the entire T — ji plane. For /x = 0, 
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we then have 



7r 



P^ = 3 — T^ 
90 



for the pressure of an ideal massless pion gas, to be compared to 



P,^37^T^-P. 



:i7) 



(18) 



for that of an ideal massless plasma of two quark flavours; here B denotes the bag pressure 
specifying the difference between coloured and physical vacua. Equating the two pressures 
(see Fig. H^) gives 

^^(^) = (3i^) ^ ' ^''^ 

so that by suitably tuning the bag pressure, we can obtain reasonable values for the 
transition temperature T^. 

For T = 0, the pressure of extended nucleons is given by 



PBifi,T = 0) 



l-nB(T = 0,/i)K 



(20) 



where hb denotes the corresponding density (eq. ( ITT]) ) and Pb the pressure of pointlike 
nucleons of mass M [19], 

n(/i,T = o) = 

1/2 






-7 



2U, 



H7)- 



M 
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The corresponding expression for massless quarks is 



P,(/i,T = 0) 



^(i' ~^' 



(21) 



(22) 



where we have converted the quark baryochemical potential to that of nucleons, fiq = ^/3. 
Again we equate the two pressures (see Fig. Hb) to obtain a critical baryochemical potential 
fic{B, Vc). Here as well one can try to tune bag pressure and hard core volume to obtain 
a reasonable threshold. 




(a) (b) 

Figure 4: Pressure comparison at /x = (a) and at T = (b) 



However, our point in presenting such toy model hadron-quark pressure comparison is 
not a quantitative determination of Tc{B) and fJ.c{B, Vc). We only want to illustrate that 
here in the low baryon density region, the confining bag pressure (as counterpart of the 
hadron size in percolation) provides the transition and sets the transition scale, whereas 
at low temperature and high baryon density, this is achieved by two competing effects, 
the confining bag pressure and the hard-core repulsion. 

Obviously, these consideration are a gross oversimplification, since we have neglected the 
resonance interactions of both pions and nucleons, as well as the role of strange mesons 
and baryons. We shall include these in the following. 

4 The Hadronic Resonance Gas 

For vanishing or low baryon number density, the interactions in a hot hadronic medium 
are resonance dominated, and hence the system can be described as an ideal gas of all 
possible resonance species [20,21], contained in an overall spatial volume V. The grand 
canonical partition function for such a gas is given by [22, 23] 

In Z(T,/i,/i5,^)=ln ZM{T,V,fis) + In Zb{T, fi, fis,V) (23) 

where the first term gives the meson and the second the baryon contributions to the par- 
tition function. Baryon number and strangeness are accounted for by the corresponding 
chemical potentials fi and /x^, respectively. For the meson contributions, we sum over all 
possible states i 

In ZMiT,V,fxs)= E In ^l,(T,V^,/i5), (24) 

mesons i 

with 

VT °° 
In Z^(T, V, fis) = d,—^ m^ ^ n~'K2{nm,/T) e^^^iMs/T^ (25) 

^^ n=l 

where di specifies the spin-isospin degeneracy and Si the strangeness of the state i. For 
the baryon contributions we have 

\n ZB{T,^i,fis,V)= E \n Z's{T,fi,fis,V), (26) 

baryons i 



with 



VT °° 

In Zi,(T,/i, V) = d^ — m\ E(-l)"^' ^"' i^2(nm,/T) e-(B'^+Si/^s)/T^ (27) 

^^ n=\ 

where Bi is the corresponding baryon number of the state i. The above form incorporates 
Bose-Einstein and Fermi-Dirac statistics; the first term of the sums in Eqs. fl25l) and (p71) 
is the Boltzmann limit. In all actual calculations we enforce both baryon number and 
strangeness conservation, with a vanishing overall strangeness. 

The resonances are confined (colour singlet) qq or qqq states of hadronic spatial size. As 
above, we consider the transition from a confined to a deconfined medium to occur when 
the hadrons as little bags fuse into one big bag, the quark-gluon plasma [24]; the critical 
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transition density is given by eq. ([8]). If we ignore for the moment the contribution 
of baryons and antibaryons, we can estimate the transition curve in the resonant, low 
baryon-density region by hmiting the meson density hm-, obtained from eq. fl25p . 



nM{T,fi) 



1.22 



0.57 fm 



-3 



(2^ 



This equation can be solved using a resonance gas summation code, summing over all 
meson states up to mass 2.5 GeV [22]; it leads for yU = to a transition temperature 
Tc ~ 177 MeV, a value considerably lower than that of pure pion gas case in the previous 
section 



Just as in the toy model, however, the temperature obtained from eq. (1251) is presumably 
somewhat too high, since we had required the box to be filled by mesons alone, while in 
fact baryons and antibaryons can contribute to forming a connecting cluster of hadronic 
matter. But including the density of baryons and antibaryons in the limiting relation (!28l) 
gives at /i = a limiting temperature Tc ~ 170 MeV, only slightly lower, so that the role 
of the baryons and antibaryons in establishing the boundary is again quite small. 

If we neglect the interrelation of strangeness and baryon number due to associated pro- 
duction, the mesonic part of the partition function is independent of fi, and hence eq. 
( 1251) leads to a constant T^. for all /x. Taking associated production into account, however, 
implies with increasing /x an increasing density of strange mesons and thus a (slightly) 
decreasing temperature. The resulting behaviour is shown in Fig. [51 
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Figure 5: Transition line in the resonance regime 

The deconfinement transition in the resonant region is conceptually a consequence of 
hadron size and clustering, leading to the fusion into one connected volume [24], and it 
occurs, as we saw above, already for a pion gas, though at a higher temperature. Thus 
the existence of the limit is due to the basic hadronic size [25], the actual value of the 
transition temperature to the scale specifying the resonance spectrum [21]. 

An interacting hadron medium can be replaced by an ideal resonance gas only if the 
interactions are dominantly of resonance nature. As the baryon density increases, however, 
there are more and more non-resonant contributions, and in the limit of low temperature 
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T and large baryochemical potential /i, i.e., for cold and dense nuclear matter, the non- 
resonant nuclear forces are the dominant interaction. 

We therefore assume as above that in the region of large baryon density, the jamming 
of baryons with a hard core defines the limit, so that the relevant relation is eq. (TT2|) . 
To connect this geometric argument with thermodynamics, we again follow the van der 
Waals approach (ITT!) for baryons and antibaryons, with Ve = 2Vc and Re = 0.4 fm. The 
resulting hard core transition curve is shown in Fig. El together with the bag fusion curve. 
The T — fi plane thus shows two distinct regimes: at low /i, hadron percolation through 
bag fusion, at large /i a first-order mobility or jamming transition. 
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Figure 6: Bag fusion vs. hard core transition lines 

As already mentioned above, the separate calculation of a bag fusion and a hard core curve 
is only an attempt to arrive at a schematic picture as long as we don't have continuum 
percolation studies for a (/z-determined) mixture of permeable spheres with spheres having 
a hard core. The curve expected from such a calculation will presumably lead to a 
transition from one regime to the other at somewhat lower values of fi. In fact, the 
baryon density becomes equal to the meson density at yU ~ 0.45 GeV, and therefore this 
may well be close to the point at which the change in the nature of the transition occurs. 

5 The Nature of the Transition 

Our considerations up to now specified the limit of hadronic matter, defined as the point of 
disappearance of the vacuum as a large-scale feature. This point was determined through 
percolation studies, and the percolation limit is in general not a thermodynamic phase 
transition. Percolation can thus naturally provide a way to produce a rapid cross-over 
not associated with any singularity of the partition function [26,27]. It should be noted, 
however, that for spin systems, thermal critical behaviour can be formulated in terms 
of percolation [28-30]. It seems possible to extend this to gauge systems, and first such 
studies relate the onset of deconfinement at /i = to Polyakov loop percolation [31], 
analogous to the onset of magnetisation as the percolation of spin clusters. Here it is 



10 



the onset of large scale disorder, i.e., of the vacuum, which induces critical behaviour. 
From the confined side, we thus have hadronic bag fusion leading to the disappearence of 
the physical vacuum, while on the deconfined side, formation of disordered clusters in an 
ordered medium correspond to the appearence of the vacuum. Based on the spin-gauge 
universality [32], deconfinement as Polyakov loop percolation could occur as first (SU(3)) 
or second order (SU(2)) phase transition, corresponding to the spontaneous breaking of 
a global center Z3 or Z^ symmetry. 

In the case of hard core percolation, a connection to thermodynamic critical behaviour 
has also been discussed [15]. If a system with hard core repulsion between its constituents 
is in addition subjected to a density- dependent negative background potential, first order 
critical behaviour can appear. A classical case is the van der Waals equation. The pressure 
in our hard core medium is given by 

1 -n[T,fi)Ve 

where Pq denotes the pressure of pointlike constituents. The density n of the hard core 
constituents is given by eq. (fTTI) . and Ve again denotes the random dense packing volume. 
If we add to this purely repulsive form a density dependent attractive term 

Pvdw(7',/u) = P{T,n)-an'^, (30) 

with constant a, we obtain the van der Waals equation of state with a first order phase 
transition ending in a second order critical point specified by the parameters a and Vg. 

Let us now recall how the transition is treated in the Nambu-Jona-Lasinio model (NJL). 
Here the thermodynamic potential Q{m) depends on the dynamically generated mass 
m ~ (ipip), which serves as an order parameter for chiral symmetry. By analyzing the 
minimum conditions, 

dQ d^n , , 

dm o^m 

at fixed T and yU, one studies the onset of chiral symmetry breaking and the order of 
the transition. At T = and /i = 0, the mass is generated through the attractive self- 
energy interaction ~ GIi/jiIj)'^, which results in chiral symmetry breaking. In a thermal 
medium of vanishing baryon density, an increase in the temperature produces fluctuations 
which eventually overcome the attraction and thus lead to the chiral symmetry restoration 
transition. At large baryon density and low temperature, the attraction has to compete 
with the repulsion due to the Fermi statistics of the quarks. As above, this competition 
brings about a first order phase transition, which persists in some temperature range 
< T <Tc] here Tc is essentially the point at which the Fermi repulsion is fully overcome 
by the thermal fluctuations. The mechanism for the large /i behaviour in the NJL model 
is thus conceptually very similar to a van der Waals pattern [33]. Note that in addition 
to the mentioned mechanisms, at sufficiently low T there will presumably be qq pairing 
into bosonic diquarks, leading to condensation and colour superconductivity. We do not 
consider this here, although it can certainly modify the resulting phase structure [34] . 

Since the NJL model does not contain gluons, it also does not include any confinement- 
deconfinement transition. The latter is implemented by introducing an interaction term 
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with Polyakov loops [12,35]. The resulting description leads in the T — fi plane to a 
deconfinement curve determined by the Polyakov loop contribution, which at some (large) 
/i intercepts the chiral curve [11,35], in a form very similar to that of Fig. O 

6 Summary 

We have argued that as function of T and fi, hadronic matter finds itself in two distinct 
regimes. At low baryon density, the behaviour of the system is governed by resonance 
formation and clustering, with hadronic size and resonance spectrum as relevant param- 
eters. The interaction here is essentially attractive. At high baryon density, in addition 
to this, there is a repulsive contribution, on the confined side as nuclear repulsion, on 
the deconfined side as Fermi repulsion between quarks. In our approach, the resulting 
limit of hadronic matter is in the low baryon density region determined by the percola- 
tion of permeable (overlapping) mesons and at high baryon density by the percolation of 
hard-core baryons. In the latter case, the competition between repulsion and clustering 
can provide a first order phase transition. In mesonic percolation there is only clustering; 
while in general not resulting in thermodynamics critical behaviour, it can in specific 
cases (depending on details of the dynamics) also result in first or second order phase 
transitions. 
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